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ON  THE  INTERACTION  OP  SHOCK  WAVES  IN  WATER-SATURATED 
EARTH  AND  IN  WATER 


G.  M.  Lyakhov 
(Moscow) 


On  the  basis  of  an  approximation  of  the  dependence  p(V)  by 
a  piecewise  linear  function^^we  gi-ve  an  approximate  solution  to  the 
problem  of  the  interaction  of  a  two-dimensional  nonstationary 
shock  wave  with  the  boundary  of  media,  a  shifted  boundary,  and 
with  a  boundary  of  a  region  of  given  pressure. 


1.  The  Interaction  of  a  shock  wave  with  a  boundary  of  media 
or  a  shifted  boundary.  In  water,  petroleum,  water  with  air  bubbles, 
water- saturated  ground,  and  certain  othei*  media  during  shock 
conq>resslon ,  as  is  known  [1-4],  the  pressure  may  be  viewed  only  as 
a  function  of  the  specific  volume  V 

r^p{y>.  (1.1) 


Pig.  1. 
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The  relationship  p(V)  is  nonlinear,  but  it  is  identical  with 
a  rise  oi  fall  in  pressure.  Thus  p'  <0,  and  p"  >  0.  In  such 
media,  called  nonlinearly  elastic,  the  basic  equations  of  gas 
dynamics,  expressing  the  laws  of  conservation  of  mass  and  quantity 
of  motion,  cind  also  the  relationships  on  the  shock  wave  front, 
conforming  to  these  same  laws,  together  with  (1.1)  form  a  closed 
system  and  allow  us  to  view  the  propagation  and  Interaction  of 
shock  waves  without  using  a  thermodynamic  relationship.  A 
thermodynaimic  relationship  remains  in  force.  Using  it,  we  may 
determine  the  entropy  losses  on  the  wave  front. 

The  propagation  of  a  stationary  wave  an  a  medium  with  variable 
sign  p"  was  examined  by  G.  I.  Barenblatt  [5]  on  the  basis  of  [6]. 
The  propagation  of  a  nonstationary  wave  was  examined  in  [?]. 

In  the  segment  h  =  0  when  t  =  0,  let  the  pressure  from  the 
shock  increase  from  Pq  =  0  to  and  fall  according  to  a  certain 

law 

(1.2) 

At  distance  h*  from  the  Initial  segment  the  medium  borders  on 
the  second  medium.  On  the  border  there  may  be  a  barrier  made  of 
incompressible  macerlal,  the  mass  of  which  is  equal  to  m  for  a 
unit  of  area. 


In  a  system  of  LaGrange  coordinates  (mass  h,  time  t)  a  solution 
of  the  equations  of  gas  dynamics  (cf. ,  for  exsunple,  [7]) 


dp 


ST 


0, 


du  _dV 
dh  IT 


0 


in  the  case  of  a  linear  dependence  between  the  volume  and  the 
pressure  has  the  fom 

—  Al)  +  /■,  (*  +  Al) 

M-pe.  p-:i/K)  (1.3) 

u  -  Ifiih-  Al)  -  /,  (fc  +  At)] 
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Fig.  2 


Here  and  further  u  Is  the  velocity  of  particles,  A  Is  the 
acouAtlc  resistance  of  the  nedlum,  and  c  Is  the  velocity  of  sound. 
Functions  F^  and  F2  are  determined  by  the  Initial  and  boundary 
conditions • 

Let  us  approximate  the  curve  p  ■  p(V)  In  the  first  medium  by 
the  two  straight  lines  (Fig.  1): 

When  o</<^ 
when 


The  niiob'er  of  approximating  links  may  be  Increased  depending 
on  the  required  dccuracy  of  the  calculation.  Prom  segment  h  ■  0 
the  front  bf- the  shock  wave  will  spread  (Fig.  2).  The  solution  In 
area  1  (elastic  wave)  we  will  obtain  In  the  form 

(1.5) 

The  form  of  function  f(t)  here  and  henceforth  Is  determined  by 
equation  (1.2). 

At  moment  t  ■  h*/Aj^  the  front  of  the  wave  reaches  the  media 
bbOndarles.  Let  In  the  second  medium  the  acoustic  resistance  A* 
lAd  ^  Us  ^  In  the  reflection  of  the  wave  the 
j>ressure  Increasllj  Ifid  \i\x9  acoustic  resistance  of  the  first  medium 
becomes  equal  to  Il2«  When  t  •  h*/Aj^  the  velocity  of  the  front  of 
the  reflected  wave  h*  ■  -A2. 
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the  ohange  In  pressure.  For  example,  in  water  with  the  fall  in 

2 

pressure  on  the  front  from  250  to  1  kg/em  the  velocity  of  the 

front  D  according  to  data  [2]  Is  reduced  from  1530  to  1500  m/s. 

In  a  water-saturated  ground  [4],  If  before  the  wave  front  Pq  ■ 

2  2 
•  25  kg/cm  ,  a  fall  In  pressure  on  the  front  from  75  to  50  kg^/cm 

Is  aocoinpanled  by  a  fall  In  the  velocity  of  the  front  by  0;08  of 

the  Initial  value  of  D  with  a  content  of  entrapped  air  equal  to 

0.01  of  the  total  volume  of  ground,  and  less  than  by  0.001  D  In 

ground  apt  containing  air.  Therefore,  we  may  assume  that  h*  ■  -Ag 

In  the  entire  segment  from  h  .h*  to  h  ■  0. 

The  solution  In  region  2  (see  Fig.  1)  Is  determined  from  the';  ; 
relationship  on  the  front  of  the  reflected  wave 

and  from  the  condition  on  the  media  boundary.  The  velocity  of  the 
medium  particles  touching  the  boundary  on  both  sides  Is  equal  to 
the  velocity  of  the  boundary  ^(t). 

Taking  Into  account  the  relationship  on  the  front  of  the 
reflected  wave  and  (1.3),  we  find  the  equation  of  motion  of  the 
boundary 

(1.6) 

and  the  solution  In  region  2  In  the  form 
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(1.7) 


#(*.  0 


■"■■  - T, - ) 


The  solution  In  region  1*  Is 


k»—k  +  jCt 

— :p — 


■)• 


(1.8) 


The  form  of  function  ^(t)  Is  determined  from  equation  (1.6). 

If  on  the  boundary  of  the  media  there  Is  no  barrier,  the  velocity 
of  the  boundary  and  the  pressure  on  the  boundary  are  determined  by 
equation 


The  ratio  of  the  pressure  acting  on  the  barrier  to  the  pressure 
in  the  Incident  wave  (the  reflection  coefficient)  In  a  nonlinear ly 
elastic  BMdlua  Is  n  (Aj^  *  A2)M^  >  2;  and  in  a  linearly  elastic 
Biedlum  when  •  A2  we  obtain  n  "  2. 

If  the  barrier  Is  absent,  then  n  ■  In  a  nonllnearly  elastic 
nediun  and  n  "  >12  a  linearly  elastic  medium;  here  for  each 
medium  we  will  get 


respectively. 

In  the  absence  of  a  barrier  the  media  boundary  moves  slowly, 
and  f(t)  decreases  monotonloally.  In  the  reflected  wave  the 
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pz*e8sure  has  a  minimum  value  on  the  front,  and  a  maximum  value  when 
h  ■  h».  In  the  second  medium  maximum  pressure  Is  attained  on  the 
wave  front,  while  minimum  pressure  occurs  on  the  media  boundary. 

In  the  presence  of  a  barrier  of  finite  mass  for  a  certain  time 
after  t  ■  t*/Aj^  the  barrier  moves  quickly,  and  ♦(t)  Increases  with 
the  Increase  in  the  argument.  The  p2*es8ure  has  its  greatest 
value  on  the  front  of  the  reflected  wave,  and  its  minimum  value 
when  h  ■  h*.  This  time  Interval  is  shown  by  Fig.  3>  Figure  4 
corresponds  to  the  moment  t2%  when  ^(t)  achieves  a  maximum.  From 
here  on  the  minimum  value  of  pressure  will  be  found  at  the  point 
which  approaches  the  wave  front  with  the  passage  of  time  (Fig.  5). 
Moments  t^,  t2«  and  t^,  to  which  Figs.  3*  and  5  correspond, 
are  represented  In  Fig.  2. 


2.  The  Interaction  of  a  shock  wave  with  a  boundary  of  an  area 
of  fixed  pressure.  The  solution  In  area  2  will  be  valid  for  the 
arrival  of  a  front  of  a  reflected  wave  to  the  segment  h  >  0. 

Let  us  examine  the  case  of  Icurge  9  and  small  h*  and  let  us 
assume  that  In  the  reflected  wave  the  acoustic  resistance  remains 
equal  to  A2.  In  the  plane  h  ■  0  p  <  p^^^.  In  the  reflection  of 
the  wave  from  segment  h  ■  0  a  change  In  the  acoustic  resistance 
takes  place  and  regions  3  and  k  are  formed  (see  Fig.  2).  In 
region  3  on  the  boundary  with  region  2  the  pressure  falls  abruptly 
to  p^y.  In  region  4  a  further  decrease  In  pressure  occurs  with 
the  growth  In  time .  On  the  boundary  of  regions  2  and  3 

(2.1) 

where  F2  corresponds  to  region  2,  and  and  ^2  correspond  to  region  3. 
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Prom  here  In  region  3 


+  ■><1  .rMt  +  At){h—Ail)  .  Mi  +  ^,)U*  1 


^  _  rg(*~^tO  .  2At*  +  AiA,  +  A*  .  .1 

+  I  AT-- Ax  +  ~-KrAx{Ax-Aii  *  J 


At—  Ax 


At 


Ax  +  At 

1  Ax  ■ 

AxAt 

2i4t 


A%X% 


*•] 


(2.2) 


The  solution  In  region  3  Is 

+  ®i)  . 


In  segment  haOp>f(t),  and  on  the  boundary  of  regions  3 
and  M  the  condition  Pjj  -  Aj^Ujj  ■  p^  -  A^^u^  ■  2*2  Is  fulfilled.  Prom 
here  It  follows  that  In  region  M  the  function  Is  the  same  as 
In  region  3i  and 


/Axt-h\ 

r^i*  —  Ax+  At  . ,1 

V  Ax,  j 

1  ±At  Aaias  SUt  ^ 

L  A  AxAt  *  J 

+ ^•»L'iir+:ir  ’unsr^T  • 


(2.3) 


The  solution  In  region  4  Is 

f.Oi  +  ct,,  (2.4) 

The  shock  wave  Is  reflected  from  the  boundary  of  the  area  of 
fixed  pressure  In  the  form  of  a  vacuum  wave. 

3.  The  Interaction  of  a  vacuum  wave  with  a  media  boundary, 
the  moment  when  the  vacuum  wave  front  reaches  the  media  boundary 
I*#il.0n8  5  and  Zf  fpe  *1  and  determine  the  flow  In 

|td|idn  Si  and  let  |hi|  dttermlne  the  flow  In  region  2*.  Prom 
thp  condition  on  bpl^dai^les  !•»  2*,  and  3»  5  we  obtain  Fg  ■  0, 
while.  8^  Is  the  same  function  as  in  region  3*  In  region  5  we 
have  pg  +  Aj^Uj  ■  region  2»  we  have  pg*  -  Aeug*  ■  0. 
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The  equation  of  motion  of  the  barrier  is 


"ti  (t)  +  Ml  +  r  W  -  M»i  (**  -  /<i«)  -  0 


(3.1)^ 


The  pressure  against  the  barrier  from  the  side  of  region  5  Is 


the  pressure  from  the  side  of  region  2  is 


(3.2) 

(3.3) 


If  the  barrier  Is  fixed,  and  when  h  ■  0  the  pressure  Is  fixed 
In  the  form 


(3.4) 


then  the  pressure  against  the  barrier  In  region  5  Is 


Ml  «f  ^>)*  Pggfg  /ZAt  +  Ai 


(3,5) 


At  the  moment  ■  h»(2A2  +  k^)/k^k2  **®8lon  7  Is  formed.  In 
this  region  the  function  Is  the  same  as  In  region  4.  If  in  the 
segment  h  ■  0  the  pressure  is  fixed  In  the  form  (3.4),  then  when 
h  ■  h®  In  region  7 


Ol  (V  -  ^,1)  «  -  ((/»,  Ml  +  + 


T7r» 


ATAiiAr+A^- 


(3.6) 


If  the  mass  of  the  barrier  Is  finite  here,  then  when  t  ■  t„ 

the  velocity  of  the  barrier  may  not  change  abruptly,  arid  therefot'e 

It  Is  determined  by  equation  (3.3).  By  using  (3.6)  and  taking 

into  account  that  ^(h*/A)  ■  0,  we  will  obtain  the  pressure  against 

the  barrier  with  t  ■  t„ 

m 

p  -t  an,  (A*  Ait^ —Alt  (r*,)  <«2""  ^  4t|  /(,  u  ^ oi**  (*"»)  ^  ®  .  (3.7) 
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The  pressure  is  negative.  The  medium  attains  a  velocity 
directed  from  the  barrier  to  the  side  of  the  segment  h  -  0,  and 
breaks  away  from  the  barrier.  Further  on  the  pressure  from  the 
side  of  the  second  medium  causes  a  reduction  In  the  velocity  of  the 
barrier  and  the  formation  of  contact  with  the  first  medium.  The 
pressure  arising  here  Is  less  than  In  the  examined  regions. 

The  author  expresses  his  thanks  to  S.  S.  Grigoryan  and 
N.  I.  Polyakova  for  their  discussion  of  the  article. 
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By  using  a  piece-wise  linear  approximation  for  the  volume  as 
a  function  of  the  pressure,  approximate  solutions  are  given 
for  three  problems.  The  first  problem  is  that  of  the  inter¬ 
action  of  a  shock  wave  with  the  boundary  of  the  medium  or  with 
a  moving  barrier.  The  propagation  of  a  stationary  wave  in  a 
region  where  the  second  derivative  of  the  pressure  with  respect 
to  the  volume  is  variable  in  sign  has  been  discussed  by  0.  I. 
Barenblatt.  Propagation  of  a  non-stationary  wave  was  considered 
in  an  earlier  paper  of  the  present  author.  At  a  given  section 
there  is  a  Jump  discontinuity  in  the  pressure  and  at  the  section 
given  the  first  siedium  is  bounded  by  a  second.  On  this  boundary 
there  can  be  a  barrier  of  an  incompressible  material.  Choosing 
h  and  the  t  as  Lagrangian  coordinates  the  equations  of  gas  dynam¬ 
ics  can  be  solved  generally  and  the  solution  applied  to  various 
regions  in  which  the  problems  of  the  paper  arise,  m  the  first 
pr^lem  we  are  concerned  with  the  propagation  of  a  wave  from 
the  given  section,  its  reflection  from  the  section  given  and 
the  propagation  into  the  second  medium,  wnlch  depends  on  the 
velocity  of  the  barrier  between  the  two  media.  Expressions  for 
the  pressure  and  velocity  are  given  in  each  case.  The  second 
problem  is  the  interaction  of  a  shock  wave  with  the  boundeury  of 
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